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We investigate the combined effect of the repuisive vector interaction and the positive eiectric 
chemical potential on the chiral phase transition by considering neutral color superconductivity. 
The chiral condensate, diquark condensate and quark number densities are solved in both two- 
flavor and two-plus-one-flavor Nambu-Jona-Lasinio models with the so called Kobayashi-Maskawa-'t 
Hooft term under the charge-neutrality constraint. We demonstrate that multiple chiral critical- 
point structures always exist in the Nambu-Jona-Lasinio model within the self-consistent mean-field 
approximation, and that the number of chiral critical points can vary from zero to four, which is 
dependent on the magnitudes of vector interaction and the diquark coupling. The difference between 
the dynamical chemical potentials induced by vector interaction for u and d quarks can effectively 
reduce the Fermi sphere disparity between the two flavors of diquark pairing. Thus the vector 
interaction works to significantly suppress the unstable region associated with chromomagnetic 
instability in the phase of neutral asymmetric homogeneous color superconductivity. 
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I. INTRODUCTION 



It is generally believed that QCD exhibits a rich phase structure in an extreme environment such as at high 
temperature and high baryon chemical potential. For the chiral phase transition, it is a widely accepted view that the 
critical point(s) (i.e. the end point of the first-order phase boundary) should exist at finite temperature and density. 
Recently, a schematic T-\i phase diagram with one critical point was extensively adopted in the literature [J, [U, 

In the last decade the color-superconducting (CSC) phase has attracted a lot of theoretical interest and prompted 
extensive studies of dense and cold quark matter [1, S B 13- At asymptotically high density that justifies the 
perturbative QCD calculations, the color-flavor locked (CFL) phase Q has been established as the ground state of 

■ quark matter. For the intermediate density region which may exist in the core of a compact star, the nonperturbative 
features of QCD play a more important role in the phase structure of QCD, and both CFL and non-CFL CSC phases 

■ may appear in this region. 
The possible emergence of a CFL or a non-CFL CSC phase around the chiral transition boundary at low tem- 
peratures should affect the chiral phase transition, and hence the interplay between the chiral condensate and the 
diquark condensate may result in an unexpected phase structure of QCD. Such an example was first presented in 'SH\ 
in a two-flavor Nambu-Jona-Lasinio (NJL) model with the vector interaction included: It was found that the repul- 
sive vector interaction can lead to a two-critical-point structure in the T-fj, phase diagram of QCD, which is due to 
the fact that the density-density correlation induced by the vector interaction effectively enhances the competition 

1 between the chiral and diquark correlations while weakening the first-order chiral restoration [J 

E3, EL We note 

that the renormalization-group analysis 13[ , the chiral instanton_anti-instanton molecule model [l4| , and the 
truncated Dyson-Schwinger model of QCD |15j all support the existence of the vector- vector four-quark interaction. 
The vector-vector interaction may be responsible for the vacuum properties of vector mesons in low-energy effective 
theories of QCD[H [U El]. 

In a quite different context, a realization of a similar two-critical-point structure of the QCD phase diagram has 
been recently conjectured [19j |: That is, the UA(l)-breaking vertex may induce a new critical point. This conjecture 
is based on a general Ginzburg-Landau theory constrained by QCD symmetries. There, the three-flavor anomaly 
term generates the cubic coupling between the chiral and diquark condensates. It has been argued that the resultant 
crossover of chiral restoration at small temperatures embodies the hadron-quark continuity hypothesis [20[. Further 
investigation is necessary to clarify whether the new critical point really exists or not with reasonable parameter. 

Beside the above two cases, another mechanism which can induce multiple critical-point structure was demonstrated 
in [2l| . where it was first disclosed that the positive electric chemical potential fi e required by the charge- neutrality 
constraint plays a similar role as the repulsive vector interaction on the chiral phase transition in a four-quark 
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interaction model: In a simple two-flavor NJL model, the same two-critical-point structure found in Q was obtained. 
In addition, positive \x e also represents the magnitude of disparity between Fermi spheres of u and d quarks when 
taking into account the local charge-neutrality constraint. For an asymmetric homogeneous system including two 
different flavor quarks with the mismatched Fermi spheres, the energy gap of the Cooper pairing between these two- 
flavor quarks can increase with temperature. This unusual behavior of the energy gap is due to the smearing of the 
Fermi surface by temperature. For the two-flavor neutral CSC phase, this unconventional thermal behavior of the 
diquark condensate can lead to a special competition between the chiral condensate and the diquark condensate, which 
may be enhanced with increasing temperature. For some model parameter regions, this abnormal competition induced 
by fi e can result in a nontrivial three-critical-point phase structure plj . The competition of two order parameters 
with the external constraint (s) should also suggest a general mechanism for a realization of a multiple critical-point 
structure, which may have some implications to study the phase transitions in condensed matter physics. 

Note, however, that a very simple two-flavor NJL model with only scalar quark-antiquark and diquark interaction 
was used in Ref. [2l|: If the model produces a relatively large vacuum quark mass, no multiple critical points appear 
in such formalism even though the charge-neutral constraint is taken into account for the CSC phase. For a more 
realistic situation, both the vector interaction and the strange quark degree of freedom should be taken into account. 

Because in the chiral phase transition similar roles are played by the repulsive vector interaction and the electric 
chemical potential under the neutrality constraint, one can expect that the chiral restoration will be weakened more 
significantly when simultaneously taking into account both the ingredients. In this case, it is possible that the multiple 
critical point structures shown in [2l| will become more robust, which we will show in the present work. In addition, 
for the neutral quark matter system, because of the quark density discrepancy, the vector interaction should also give 
different contributions to the dynamical chemical potentials for u and d quarks. This may have important influence on 
the property of the asymmetric homogeneous CSC phase which suffers from the so called chromomagnetic instability 
characterized by the negative Meissner mass squared [22J. This is a significant observation and will be shown to be 
the case in the following sections. 

There have been extensive studies within the NJL model on the phase diagram of dense and locally neutral three- 
flavor quark matter with strange quarks explicitly included [23l. [2^|: In the strong coupling case, the chiral breaking 
phase at low temperature is bordered by the 2CSC phase, while for the weak coupling case it is surrounded by both 
normal quark matter and a gapless phase called g2CSC. However, no new critical point was found in [231 [24] , in 
contrast to the result reported in [21|. Notice, however, that the vector interaction was not included in [23l. |24| . 
and hence only a strong first-order chiral restoration happened in the low temperature region. We can expect that 
the qualitative feature of the phase diagram, even for the three- flavor quark matter reported in [23|, [24[ , may change 
when the repulsive vector interaction is taken into account together with the charge-neutrality constraint. Again, we 
shall see that this is actually the case, and there appear multiple critical points in the phase diagram of three-flavor 
quark matter. 

The paper is organized as follows. In the next section, a nonlocal two-flavor NJL model is used to investigate 
the (T, /j,) phase structure of QCD, with both vector interaction and charge-neutrality being taken into account. The 
extension of the work to the two-plus-one-flavor case with the so-called Kobayashi-Maskawa-'t Hooft(KMT) term(25jis 
presented in Sec. III. The final section is devoted to a summary and concluding remarks. 



II. TWO-FLAVOR CASE 



We will demonstrate the joint effect of the vector interaction and neutral electric-charge constraint on the QCD 
phase diagram within a nonlocal two-flavor NJL model. 



A. Model 



NJL-type models have been extensively used to investigate the CSC phase transition at moderate and large den- 
sities @: See the previous study without the CSC P, [ljj, IH, HtJ- The advantage of a NJL model is that it can 
investigate the interplay between the chiral condensate, the di quar k condensate and the quark number density on the 
same footing. For the two-flavor case, a nonlocal NJL model [28|, [29|, [30, [3l|, [32|, [33|, [34[ is adopted here, which has 
the Lagrangian density, 

3 

- G v [{q(x)Y L T iq (x)) 2 + (q(x)t^ l5 r iq (x)) 2 

i=0 

+ G D ^2 \q{x)ibT2^Aqc{x)] [qc(xh5T 2 \Aq(x)} , (l) 

A 



C = V m)ip + G s (q(x)q(x)) + (q(x)ij 5 Tq{x)Y 



3 



Parameter set 


m(MeV) 


GgA 2 


A (MeV) 


-{uu} 1/3 (MeV) 


M(p = 0) (MeV) 


Set 1 


5.01137 


2.64310 


600.271 


248.440 


400 


Set 2 


4.92671 


2.51088 


617.968 


249.906 


367.5 


Set 3 


4.70805 


2.35908 


649.168 


253.699 


330 



TABLE I: Model parameter sets for two-flavor nonlocal NJL. 



where 

q(x) = J dy 4 f(x - y)i>(y) , q c (x) = J dy A f{x - y)tpc{y) , ipc = G^ T (2) 

and C — «7o72 stands for the Dirac charge conjugation matrix. The three coupling constants, namely Gs, Gy, and 
Gd, belong to the scalar mesonic channel, the vector mesonic channel, and the diquark channel, respectively. The 
current quark mass matrix is given by fh = diag(m u , md) in two flavors, and we shall work in the isospin symmetric 
limit with m u = m<j = rn. We note that A^'s are the antisymmetric Gell-Mann matrices (i.e. A runs over 2, 5, 7 only) 
for the color SU(3) group, while tq and the f's are the unit matrix and Pauli matrices in flavor space, respectively. 
In contrast to the scalar interaction, the vector part has U(2)xU(2) flavor symmetry, and hence the vector terms can 
be decomposed into ~ (uj^u) 2 + (dj^d) 2 without the flavor mixing term like u^^udj^d, although there are terms 
like u'y^dd'jfiU. 

In contrast to the local NJL model, Lagrangian (JTJ) is formulated with a nonlocal interaction, which is controlled by 
a form factor f(x). This type of model can be considered as a special case of the truncated Dyson-Schwinger equation 
with a separable effective gluon propagator or an instantaneous nonlocal chiral quark model. The main purpose of 
adopting such a nonlocal interaction model in our study is to deal with the one-loop ultraviolet divergence, especially 
for the calculation of the Meissner mass squared at finite temperature [35]. For convenience, we follow the model 
parametrization in [34| and use the so-called Lorentzian form factor with the form 

/ 2 (p = lPl)=g(p) = 1 + ( 1 f)2o , (3) 

where f(p) is the Fourier transformation of the form factor fix). In Eq.([3|), a is a dimensionless parameter and A 
stands for the scale parameter of the model. To check the sensitivity of the main results on the model parameter 
choice, three sets of model parameters are adopted in this paper, which are listed in Tabel.)^) 1 By choosing a = 10 
2 , the other three model parameters, namely A, Gs and m, are determined by the vacuum physical quantities of the 
pion mass M v = 135MeV, the pion decay constant f n = 92.4MeV, and the quark condensate — (uit) 1 / 3 w 250MeV. 
Note that the value of the current quark mass is around 5MeV and the Gell-Mann-Oakes-Renner relation holds well 
for these parametrizations of the model parameters (34|. 

In the nonlocal NJL model, the produced constituent quark mass is momentum dependent. Table ((XJ) shows that 
the vacuum constituent quark mass M(p = 0) ranges from 330 MeV to 400 MeV for the listed three sets of model 
parameters which almost reproduce the same vacuum physical quantities. Although the vacuum constituent quark 
mass is not an observable, all the values of M(p = 0) in Table (P) are phenomenologically acceptable because the mass 
3M(p = 0) is larger than the nucleon mass in the vacuum, which is a bound state of three quarks. 

Because there are no reliable constraints on Gy and Gd within the two-flavor NJL model, these two model pa- 
rameters are taken as free parameters in our treatment. The standard ratios of Gy/Gs and Gd/Gs from Fierz 
transformation based on the local color current-current interaction are 0.5 and 0.75, respectively. In addition, in the 
molecular instanton liquid model, the ratio Gy /Gs is about 0.25. Combining this point with the Fierz transformation, 
the reasonable value of Gy /Gs may be located in the range from to 0.5. In the following numerical calculations, we 
will focus on the influence of the vector interaction on the phase diagram by varying Gy /Gs while fixing Gd/Gs as 
its standard value. In the literature, such a choice of diquark interaction is usually called the intermediate coupling. 



1 We mention that our results obtained in the Lorentzian nonlocal cutoff scheme are equally obtained even with the traditional sharp 
cutoff scheme. In this sense, our results, especially the realization of the multiple critical-point structures of the phase diagram are 
robust and not artificial. 

2 We have checked the main conclusion of this paper is insensitive to a or the form of g(p) if it can give reasonable vacuum properties of 
QCD. 
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B. Thermodynamic Potential with Neutrality Condition 

The grand partition function is given by 

Z = e- ny / T = [ D^e^ dx ^ c+ ^'^\ (4) 



where is the thermodynamic potential density and fi is the quark chemical-potential matrix. In general, the quark 
chemical-potential matrix fi takes the form (36j 

p, = p, - p e Q + h 3 t 3 + p 8 T 8 , (5) 

where p, is the quark chemical potential (i.e. one third of the baryon chemical potential), p e is the chemical potential 
associated with the (negative) electric-charge, and /z 3 and ps represent the color chemical potentials corresponding 
to the Cartan subalgebra in color SU(3) space. The explicit form of the electric charge matrix is Q = diag(|,— |) 
in flavor space, and the color charge matrices are T 3 = diag(|,— 1,0) and T 8 = diag(|, |, — |) in color space. The 
chemical potentials for different quarks are listed below: 

211 211 

Pru = p — 3 Me + + 3P8 , P-gu = P — ^Pa _ 2^3 3^8 j 

H rd = p + i/i e + §/i 3 + lfl S , p gd =/-!+ \He - 5M3 + \p& , (6) 

2 2 ,12 

Pbu = M — 3 Me — 3 M8 ! MM = A* + 3 A*e ~ 3 A*8 ■ 

At finite temperature and density, the Lorentz invariance is broken. The three types of four-quark interactions 
in Eq. ([T]) could develop three different condensates in a homogeneous phase: a a = (only sum over a), 

A = ((9c)a*75e a £ab39^), an d Pa = {<lal la) ( 0Iu y sum over a )i where the indies a and a refer to flavor and color, 
respectively. The last vector-type condensate appears due to the breaking of Lorentz invariance at finite temperature 
and density, which corresponds to the quark number density for flavor a. Here, we follow the common treatment for 
the two-flavor CSC phase where the blue quarks do not take part in the Cooper pairing. 

In the mean-field approximation which we will adopt, it is convenient to introduce the following two gaps and an 
effective chemical potential corresponding to the above condensates; that is, the dynamical quark mass 

M(p) = m - 2G s (a u + a d )g(p) , (7) 

the gap for CSC phase 

A(p) = A 5 ( P ), (8) 

and the dynamical quark chemical potential 

Aa(p) = - 4G v p a g(p)- (9) 

Note that the nonlocal interactions lead to all the gaps being momentum dependent, which is described by g(p) in the 
separable model. The advantage of the choice of the separable interaction is that the gap equations can be obtained 
by the variational method. We also note that the induced chemical potentials, —AGvp a g(p), for u and d quarks are 
different from each other due to the constraint of electric charge neutrality (fid > A*« an d hence pd > p u )', notice also, 
however, that they are dependent only on the respective density p u ,d, and hence the dynamical chemical potentials 
Pu,d(j>) tend to come closer because pd > p u with the common coupling constant Gy . 

Using the standard bosonization technique, the mean-field thermodynamic potential in the NJL model with the 
electron contribution takes the following form: 

n = - T H [ tS^ 1 - SmU ^ p1 + n c , (io) 



where 



Q L = 2G s (al + a 2 d ) - 2G v {p 2 u + Pd ) + - JL L* + 2^T 2 p\ + ^T 4 ) , (11) 

and the sum runs over the Matsubara frequency oj n = (2n + \)ttT and Tr is taken over color, flavor, and Dirac indices. 
The last term in (fTTj) is the contribution of the free electron gas. Note that a UV divergent counter part fic is 
introduced in ([TO]) due to the one-loop integration. 
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The inverse quark propagator matrix in the Nambu-Gor'kov formalism is given by 
with 

[Gtr 1 = 7o(^„ ± Hp)) -1-P- M(p) . (13) 
Taking the Matsubara sum, we can express the thermodynamic potential as usual as 



n(/ie, Ms, Ms, <r, i/, A;p,T) = Q L -f^J -0^{(E t - Ef) + 2Tln(l 



e 



-Ei/T 



)}, (14) 



with the dispersion relations for six quasi-particles (that is, 2 flavors x 3 colors; the spin degeneracy is already taken 
into account in Eq. ([TljO and six quasi-anti-particles. In Eq. (fTl|) , the counter part corresponds to 

d3 P 



, ( 2 -) 

where E® = Ei(M = m, A = 0, p = 0) . The unpaired blue quarks have the following four energy dispersion relations, 
Ebu = E — jlb u , Eb u = E + jlbu Ebd = E — fibd , Ebd = E + jlbd (16) 



with E = yjp 1 + M 2 (p) . In the rd-gu quark sector with pairing we can find the four dispersion relations, 

Efd-gu ~ Ea i \ {jird ~ figu) = E A ± 5fl , 
E rd-gu = EA ± \{JXrd ~ frgu) = E A ±8fl 1 

and the ru-gd sector has another four as follows: 

E ru-gd =EA± |(Ar« ~ figd) =£aT^, 
Efu-gd = Ea ± \ (&ru ~ P-gd) =£aT^, 



(17) 



(18) 



where E A = y/(E - jj) 2 + A{p) 2 and E A = y/(E + jl) 2 + A(p) 2 . The average chemical potential is defined by 

M = ^— — = ^— — = M - -g- - 2G y (p u + Pd)ff(p) + y , (19) 

and the effective mismatch between the chemical potentials of the u quark and the d quark takes the form 

5p=\{p e -AGv{ P d- Pu)g{p))- (20) 

We notice that the vector interaction coupled to the difference of the u and d quark densities tends to diminish the 
mismatch in the chemical potentials. This effect will play an important role for avoiding the color magnetic instability: 
see below for further discussions. 

For the two-flavor CSC phase, pis always remains zero due to the left unbroken color SU(2) symmetry for red and 
green quarks. Minimizing the thermodynamic potential (|14|) . we can solve the mean fields a , A, p u \d\ together with 
the chemical potential p e and ps from 

on _ on _ on _ on _ on _ on _ 

da dA 3 dp u dp d dp e dp 8 

where a — a u = ad as demonstrated in [2lj |. Since p% is very small for color neutralized two-flavor CSC matter, 
ignoring this term has a little effect on the phase structure 3 . Considering this point, Eq. (|2"5|) is then simplified as a 
series of five coupling equations. 



3 We have checked that the explicit introduction of /xg does not alter the result that the multiple critical-point structures are realized, 
although there is a slight change in the parameter window for realizing them. It is rather amazing because the value of /ig can be 
in the same order as that of the induced chemical potential near the phase boundary. The reason of the robustness of the multiple 
critical-point structures lies in the fact that such phase structures are mainly driven by the competition between the chiral condensate 
and the diquark condensate in the coexisting region while nonzero /ig has little effect on this competition. The contributions of fig to 
the red and green quark chemicals have a negative value, fig/3, while that to the blue quark has a positive value, -2fig/3 with a doubled 
absolute value, implying that the effect of fig on the quark mass and diquark condensate tends to cancel with each other. 
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V^d- 4 Pd G v9(P) 



Pu=^u- 4 Pu G v9(P) 



G v =0 
G v >0 



> P 



FIG. 1: Effect of the vector interaction on the chemical-potential disparity between the u quark and the d quark. Notice that 
the amount of the decrease from the chemical potential to the effective chemical potential jl U: d is proportional to the respective 
density p Ut d and hence the difference between the effective chemical potentials becomes significantly smaller. 



To determine the region of the unstable homogeneous CSC phase associated with magnetic instability, one need 
to calculate the Meissner masses squared which may be negative with charge neutrality. Here we adopt the same 
method as in [35} to evaluate the Meissner mass squared 

m 2 M = ^[Q(A) - n(A = 0)} B=0 , (22) 

where B has the same meaning as in [35] . In our calculation, the upper limit of the momentum integration in Eq. (|22|) 
is infinity, which helps to avoid the cut-off sensitivity for the evaluation of the Meissner masses encountered in the 
conventional local NJL model (35|. 

Note that there are three main changes induced by the nonzero vector-type quark condensates in comparison to 
the case without vector interactions [2 if ]. First, it give new negative contributions to the thermal potential, which 
favors the phase with relatively larger dynamical quark mass. This effect becomes more significant when the quark 
number density is sizable. Second, it give rise to a negative dynamical chemical potential, which can delay the chiral 
restoration towards larger chemical potential to drive the formation of the coexistence (COE) phase with both the xSB 
and the CSC phase. In the COE region, the competition between the corresponding order parameters can significantly 
weaken the first-order chiral phase transition. Third, as shown in Fig[T] and briefly mentioned above, the disparity 
between the densities of u and d quarks can effectively suppress the chemical-potential mismatch between these two 
flavors, which might partially or even totally cure the chromomagnetic instability. More details of the influences of 
the vector interaction combined with the electric-charge neutrality on the phase of QCD will be given in the next 
subsection. 



C. NUMERICAL RESULTS AND DISCUSSIONS 



In this subsection, we shall discuss the effect of vector interaction combined with the charge-neutrality and /3 
equilibrium on both the chiral phase transition and the instability of the CSC phase. Two points will be stressed 
below: In general, the model parameter window always exists in the NJL model, which favors the multiple critical- 
point structures, and the number of the critical points can be zero, one, two, three, and even four; the different 
dynamic chemical potentials induced by the vector interaction for u and d quarks can effectively shrink the unstable 
homogeneous CSC region towards lower temperatures and larger chemical potentials. 



1. Multiple critical-point structures for chiral restoration 



In the subsequent subsections we shall present numerical results and see that what is described above is actually 
the case. For convenience, we shall adopt the same notations as those in Refs. [l9|, [2l[ to distinguish the different 
regions in the T-fi phase diagram: NG, CSC, COE, and NOR refer to the hadronic (Nambu-Goldstone) phase with 
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FIG. 2: The phase diagrams for model parameter set 1 with varying Gv/Gs and fixed Gd/Gs = 0.75. NG, CSC, COE, 
and NOR refer to the hadronic (Nambu-Goldstone), color-superconducting, coexisting and normal phase, respectively. The 
boundary of the unstable region from the chromomagnetic instability is indicated by the dash-dotted curve. With the increase 
of Gv/Gs, the number of critical points changes and the instability region tends to shrink toward the lower T and higher jj, 
region in the phase diagram. 

(7^0 and A = 0, the color-superconducting phase with A/0 and a = 0, the coexisting phase with cr ^ and A ^ 0, 
and the normal phase with a = A = 0, respectively, though they have exact meanings only in the chiral limit. 

In Ref. [2lj . four types of critical-point structures were found with varying diquark coupling constant. We will 
change the value of the vector coupling Gv/Gs while fixing the diquark coupling at a standard value. The vector 
coupling constant in the vacuum may be determined by the vacuum properties of vector mesons. Since the usually 
adopted scale parameter in the two-flavor NJL model is less than the minimum vacuum mass of vector mesons, 
one could not get a reliable vector coupling constant within such formalism; on the other hand, the vector coupling 
constant should also be a function of temperature and quark chemical potential, which is also unknown at the moment. 
Therefore, to observe the possible effect of vector interaction on the phase transition, the vector coupling constant 
Gv/Gs is treated as a free parameter rather than a fixed parameter in the model at hand. 

The phase diagrams for parameter set 1 in TabeUwith different values of Gv/Gs are shown in FigO Parameter 
set 1 corresponds to a relatively large vacuum constituent quark mass, M(p — 0) = 400MeV. Figured indicates 
that chiral restoration at low temperature keeps a first-order transition, and only a single chiral critical point app ears 
in the phase diagram with vanishing vector interaction. In contrast to the case with the same Gd/Gs in [2l| . the 
first-order chiral transition is too strong to favor the appearance of the COE even though the finite fi e significantly 
moves the boundary of chiral restoration towards larger /i as shown in Figj3] With nonzero Gv/Gs = 0.32, however, 
Fig. \2)p shows that COE emerges and a three-critical-point structure for the chiral transition is realized. This result 
further verifies the conclusion in Ref. pll ] that the abnormal thermal behavior of the energy gap for the mismatched 
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FIG. 3: Chiral phase transitions for model parameter set 1: Numbers 1, 2, 3, and 4 correspond to the cases with vanishing Gv 
and no charge-neutrality constraint, vanishing Gv and a charge-neutrality constraint, Gv/Gs = 0.32 and no charge-neutrality 
constraint, and Gv/Gs = 0.32 and a charge-neutrality constraint, respectively. The standard value Gd/Gs = 0.75 is adopted 
in all cases. 

diquark pairing can give rise to two new critical points in the low temperature region. We notice that Fig. [2Jd has 
a similar phase structure to that given in Fig. 4b of plj . and the mechanism to realize this structure is understood 
much the same way: The first-order transition line E-F is the remnant of the chiral transition without the CSC phase, 
which tends to cease to exist at high temperature, while the other first-order transition line ending at G is the chiral 
transition that survives the effect of a rather strong diquark condensate at low temperatures. In addition, Fig[2j3 also 
shows that all three critical points are free from the chromomagnetic instability. 

Note that the vector interaction with Gv/Gd = 0.32 is not yet strong enough to lead to the emergence of COE 
without the help of the constraint of electric-charge-neutrality, as indicated in Fig. [31 In contrast to the case without 
both Gy and the constraint of electric charge, Fig. [3jshows that the critical chemical potential for the chiral transition 
at zero temperature is delayed towards larger fi by about 42MeV by nonzero Gv and /i e . Therefore, the combined 
influence of these two elements on the chiral transition is quite significant. 

By further increasing Gy, Fig(2Ji shows that the remnant first-order chiral transition turns into a crossover and 
the surviving first-order transition still exists with the emergence of another new end point, H. This is quite different 
from the result in Ref.pH that the surviving part of the first boundary initially vanishes with increasing Go- The 
reason for the difference is that the relatively large vector interaction has more significant impact on the location of 
the critical point E, which shifts to lower T and larger /i and eventually is eaten by the enlarged COE. This point is 
also indicated clearly in Figj3j The appearance of H is due to the stronger competition between the chiral condensate 
and the diquark condensate at lower temperatures. Note that the two critical points in FigJ5Ji are also free from 
the chromomagnetic instability. Figure^ indicates that one critical-point structure appears again in a very narrow 
parameter region of the vector coupling, where the first-order line only contains the surviving part. 

For larger vector coupling, the chiral boundary will totally become crossover, which is not shown in FigO In 
conclusion, with parameter set 1 and Gd/Gs = 0.75 , five types of critical-point structures with critical-point numbers 
1, 3, 1, 2, and are found in the model when Gv/Gs is varied in the range — 0.5. In contrast to the situation in 
Ref.[21j|, all the multiple critical points obtained with parameter set 1 are far from the unstable homogeneous CSC 
region. 

The phase diagrams for parameter set 2 with varying Gv/Gs are shown in Fig[3J Compared to parameter set 1, 
the resulting vacuum quark mass M(p = 0) is reduced and the first-order chiral transition is not very strong. For 
vanishing Gy, FigfJJi gives a similar phase diagram as Figj2^. When Gv/Gs is increased to 0.253, FigJSJa tells us 
that a new critical point, H, appears at a very low temperature. This is not surprising since the larger the diquark 
condensate in the COE phase, the more suppressed the chiral condensate in the lower T and larger \x region, where 
the phase change is a crossover. With the further increased vector coupling, FigQJ: shows that four critical points 
appear on the chiral boundary: Owing to the abnormal T dependence of the diquark condensate, two new critical 
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points denoted by F and G exist in the phase diagram. Because the vector interaction is not as strong as the case 
in FigJ2t, both the remnant and the surviving first-order chiral transition remain in the phase diagram. One can 
see from FigHb and FigHh that the critical point labeled as H is located on the border between the stable region 
and the unstable region, while other critical points are free from the chromomagnetic instability. In Figf5]i, with 
much lager Gv/Gs, the surviving first-order transition in the lower T region can not survive anymore, and only the 
remnant first-order transition in the high T region is left on the chiral boundary. Therefore, for parameter set 2 and 
the standard diquark coupling, five types of chiral critical point structures also exist, and the order of the number of 
the critical points is 1,2,4,2,0 with increasing vector interaction. 

In Table HU model parameter set 3 reproduces the smallest vacuum dynamical quark mass, and the corresponding 
first-order chiral transition is also the weakest one. Figure|5Ji shows that only finite fj, e is capable of realizing the 
COE in the phase diagram, even though the competition between the chiral and CSC correlations in the COE is 
not strong enough to lead to a multiple critical-point structure. With the help of vector interaction, the three- and 
two-critical-point structures appear in Figs. [SJd and [St, respectively. In contrast to the former two cases, only a very 
small vector interaction can lead to multiple critical-point structures. 

Since all the critical-point structures mentioned above are obtained within the range Gv/Gs < 0.35, one can 
expect that similar results may be obtained for a relatively weak diquark coupling with the vector coupling varying 
in the range < Gv/Gs < 0.5. In general, at least in the two-flavor NJL model, we can conclude that there always 
exists a region in the parameters plane of Gd~Gv which favors a multiple chiral-critical-point structure, with color 
superconductivity and electric neutrality being taken into account. 

Here it is worth emphasizing that the abnormal T dependence of the diquark condensate is a general feature in the 
COE region, even though the same behavior only happens in the CSC dominant region with a weak diquark coupling. 
The reason is that the quark Fermi surface in the COE is relatively small due to the sizable quark mass. Note that, 
the various critical-point structures demonstrated above are realized in a very simple chiral quark model of QCD with 
the physical constraint of charge neutrality. Because of the complexity of QCD, it is possible that there are more 
than one chiral critical point in the true QCD phase diagram . 

2. Suppressing the chromomagnetic instability 

It is well known that the asymmetric homogeneous g2CSC phase suffers from the chromomagnetic instability. At 
zero temperature, the calculation based on the hard-dense-loop method |22j suggests that the Meissner mass squared 
of the 8th gluon becomes negative for ^ > 1 while the 4th- 7th gluons acquire negative Meissner masses squared 
for -£ > \ j\f2. Note that without the vector interaction, the chemical-potential mismatch is just equal to /j, e /2. 
The instability of the homogeneous CSC phase should imply the existence of a yet unknown but stable phases in this 
region of the phase diagram. Such examples proposed so far include the Larkin-Ovchinnikov-Fulde-Ferrel (LOFF) 
phase and the gluonic phase [38[ . 

Because the Fermi surface is smeared by finite temperature, the homogeneous neutral two-flavor CSC phase can be 
stable in a somewhat higher temperature region [35l l39l . 1401 ] . The same thing happens for the homogeneous neutral 
CFL phase with three flavors [4lJ. It is also reported in Ref.[42] that a large quark mass and a strong coupling 
can effectively suppress the instability even at zero temperature. Here we will demonstrate that the repulsive vector 
interaction can also resolve or suppress the instability problem. The reason is very simple: According to Eq. (I20p . 
the density mismatch between u and d quarks gives a negative contribution to Sfi, which can effectively suppress the 
ratio 6fi/A. 

The ratio 5/z/A (here the ratio S/i/ A refers to 6fi(p — 0)/A(p = 0)) as a function of /j, with fixed T = 5MeV for 
different Gv/Gs is shown in FigJ5^, where model parameter set 2 is used with Gd/Gs = 0.75. One can see that 
Sfi/A significantly decreases with an increasing vector coupling. The regions of the chromomagnetic instability for 
the corresponding Gv/Gs are shown in FigO). With increasing Gv/Gs, the unstable region shrinks towards a lower- 
temperature and higher chemical-potential region. Actually, we have already seen that the same phenomenon occurs 
in the phase diagrams of Figsl^HU 

We should notice here that, as seen from Eq. (|27p . the repulsive vector interaction also suppresses the magnitude 
of the diquark condensate due to the reduced effective quark chemical potential. However, the direct effect of the 
vector interaction on S/i is more significant than that on A, in particular for finite temperature. We stress that this 
important role of the vector interaction in the CSC phase, especially for the instability problem, was first revealed 
in this work. Of course, the vector interaction may not totally remove the unstable region from the phase diagram 
unless the diquark and/or vector coupling are very large. Therefore, other mechanism may still be necessary for a 
thorough cure of the magnetic instability. 
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FIG. 4: The phase diagrams for model parameter set 2 with varying Gv/Gs and fixed Gd/Gs = 0.75. The unstable region 
with chromomagnetic instability is indicated by the dash-dotted curve. 



III. TWO-PLUS-ONE-FLAVOR CASE 



In this section, the study of the influence of the vector interaction on the chiral phase transition is extended to the 
two-plus-one-flavor NJL formalism. 



A. Model 



In this part, for simplicity, a local two-plus-one-flavor NJL model is adopted 



ocal two-plu 

was developed in the mid-1980s [HI El 0] 
interaction term and a determination term [25 
we take the same model parameters as in Ref. 
Lagrangian density is given by 



The two-plus-one-flavor NJL model 
and the most popular version includes a chiral symmetric four-quark 
in flavor space [H [46|, |47|. To compare with the previous study, 
231 by including the vector interaction channel. The corresponding 



O 

C = tW-mW + GsY, [(^V) 2 + (^7 5 A,V) 2 ] -GvY, [(^A^) 2 + (^75^) 



i=0 



+ GdY, [i>iil5e a ^e abc (^c) b p] [^cY^l^ e rsc ^] - K I det [$ (1 + 75) 1&] + det (1 - 75) V>] > , (23) 

7,c ■* 
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FIG. 5: The phase diagrams for model parameter set 3 with varying Gv/Gs and fixed Gd/Gs 
with chromomagnetic instability is indicated by the dash-dotted curve. 



0.75. The unstable region 




FIG. 6: The effect of the repulsive vector interaction on the ratio Sfi/A (left figure) and the location of the boundary between 
the stable and unstable homogeneous regions(right figure). With the increase of the ratio Gv/G s = R, S/J./A tends to become 
smaller and the unstable region in the T-/J, plane shrinks toward a lower T and higher \i region. 
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m u , d (MeV) 


m 3 (MeV) 




ifA b 


A (MeV) 


M Uld (MeV) 


5.5 


140.7 


1.835 


12.36 


602.3 


367.7 


/.(MeV) 


m^(MeV) 


m K (MeV) 


m, v (MeV) 


m^MeV) 


M s (MeV) 


92.4 


135 


497.7 


957.8 


514.8 


549.5 



TABLE II: Model parametrization of the two-plus-one-flavor NJL model. 



where the quark spinor field tp® carries color (a — r, g, b) and flavor (a = u, d, s) indices. In contrast to the two-flavor 
case, the matrix of the quark current masses is given by m = diag^(m u , nid, m s ) and the Pauli matrices in flavor space 

are replaced by the Gell-Mann matrices \i in flavor space with i — 1, ... ,8, and Ao = -\/2/3 1/. The corresponding 
parametrization of the model parameters is given in Table (|n| ). where G5, the coupling constant for the scalar meson 
channel, and K, the coupling constant responsible for the 17a (1) breaking, or the KMT term(25j. are fixed by the 
vacuum physical observables. The other two coupling constants, Gy and Gd, are still taken as free parameters in the 
following. 
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B. Thermodynamic potential for neutral color superconductivity 



In general, there exist nine possible two-quark condensates for the two-plus-one-flavor case with the Lagrangian 
(f2"3")) : three chiral condensates <j a , three diquark condensates A c , and three vector quark condensates p a , where a 
and c range from 1 to 3, which stand for three flavors and three colors, respectively. At the mean-field level, the 
thermodynamic potential for the two-plus-one-flavor NJL model including the charge-neutrality constraints, is 

3 3 3 

= n l + ^-Y,\A c \ 2 -2G v J2pl + 2G s J2al 

c—1 a — 1 a — 1 

- iK* u * d * s - ^ ^ In det ^ , (24) 

where fi; stands for the contribution from free leptons. Note that, for consistency, f2; should include the contributions 
from both electrons and muons. Since » M e and M e ~ 0, ignoring the contribution from muons has little effect 
on the phase structure. Therefore, we use fli given by the last term in Eq. pip . 

It should be stressed here that, for simplicity, the contributions of the cubic mixing terms among three different 
condensates, such as a A 2 , pA 2 , and ap 2 , are neglected in Eq. ()24[) . These terms arise from the KMT interaction which 
may or may not affect the phase structure. In particular it was argued in fl9| that the cubic mixing term between 
chiral and diquark condensates may play an important role in the chiral phase transition in the low temperature 
region. Beside the direct contribution of these cubic terms to the thermodynamic potential, the flavor mixing terms 
arising from the KMT interaction also have influence on the dispersion relationship of quasiquarks. For example, 
both the diquark condensate and the quark number density contribute to the dynamical quark mass. Therefore, it is 
a very interesting topic to investigate the possible effect of these cubic coupling terms on the phase diagram by using 
dynamic models of QCD. Leaving the discussion of this interesting problem to our future work, here we just simply 
assume that none of these mixing terms makes a qualitative difference in the phase diagram. 

Because of the large mass disparity between the strange quark and the u (d) quark , the favored phase at low 
temperature and moderate density might be 2CSC phase rather than CFL as demonstrated in the two-plus-one-flavor 
NJL model[H, [24|. These studies suggest that the strange quark mass is close to or even larger than p near the 
chiral boundary, which means that the strange quark density is considerably smaller than that of the u and d quarks. 
Therefore, for the two-plus-one-flavor case, the light quarks still play the dominant role around the chiral boundary, 
and the density mismatch between u and d quarks under the electric-charge-neutrality constraint is still similar to 
the two-flavor case. Since the main purpose of our study is to investigate the influence of the neutral CSC phase 
on the chiral phase transition by taking into account the vector interaction, we only consider the 2CSC phase in the 
following. 

The inverse quark propagator in the 2CSC phase in the two-plus-one-flavor case still takes the same form as S^ F , 
with the extended matrixes p and M in three-flavor space. The constituent quark mass is given by 

M a =m a - 4:G s <J a + 2Kapo 1 , (25) 

and the effective quark chemical potentials take the form 

2 

p u = p - AGvPu - -^Pe, (26) 

Pd = P - iGvPd + ^Me, (27) 

p s = p - 4G v p s + -p e . (28) 

The quantity p ( 5p ) still has the same form as Eq. (flT)|) [Eq. ([2H)) ] with g(p) = 1. Ignoring the mass difference between 
the u quark and the d quark (the mass difference is very small [23j), the last term in Eq. (|24|) has an analytical 
form which greatly simplifies the numerical calculation. Adopting the variational method, we get the eight nonlinear 
coupling equations 

on _ on _ on _ on _ on _ on _ on _ on _ Q , 

do u da s dA dp u dpd dp s dp e dps 

Since p$ is near zero around the chiral transition region [2^, [24[, taking it to zero makes little difference in the 
calculational result and the nonlinear equations can be reduced to 7. 
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C. Numerical calculation and discussion 

Similar to the two-flavor case, the vector interaction Gy/Gs is also taken as a free parameter in the following, and 
the diquark coupling Go /G$ is fixed to the standard value. Because of the contribution from the KMT interaction, 
the ratio GnJGs from Fierz transformation should be 0.95 rather than 0.75 when only considering the four-quark 
interaction Q- Note that in this case the effective four-quark interaction which determines the quark constituent 
mass in vacuum is Gs> = G$ — \Ka sl and the standard value of Gd/G' s should be 0.75. 

The phase diagrams with a multiple critical-point structure for different vector interactions are shown in FiglT] 
One can see that these phase diagrams are very similar to FigJH the two-flavor case with parameter set 2. This is 
reasonable since these two models almost reproduce the same constituent quark (u and d) masses and have similar 
scale parameters. Figure[7] indicates that the KMT interaction does not change the possible multiple critical-point 
structures for the chiral phase transition. 

For simplicity, the unstable regions with chromomagnetic instability are not plotted in Fig(7) The calculation of 
the Meissner mass squared in the 2CSC phase for the two-plus-one-flavor case is straightforward but complicated. 
Including the s quark should have little effect on the value of the Meissner masses calculated according to the formula 
for the two-flavor case [22j since the s quark does not take part in Cooper pairing. We can expect that the critical 
points E, F, and G should still be free from the chromomagnetic instability, as in Figs[2]and[U since the large strange 
quark may give a positive contribution rather than a negative one to the Meissner masses squared. As for critical 
point H, it may be located in the unstable region and could be safe from the instability because the relatively large 
Gv/Gs may suppress the unstable region to lower T and higher /i. 

Usually, it is argued that the instantons should be screened at large chemical potential and temperature. Therefore, 
compared to it's vacuum value, the coupling constant K is expected to be reduced around the chiral boundary. For 
smaller K, the flavor mixing effect is suppressed and the mass mismatch between the s quark and the u(d) quark 
becomes larger. Accordingly, the influence of the s quark on the chiral restoration is weakened, and the situation 
approaches the two-flavor case. On the other hand, with decreasing K, the u(d) quark mass also decreases since 
the contribution from the s quark mass is reduced. This means the first-order chiral restoration will be weakened 
when decreasing K. Correspondingly, the COE region should be more easily formed with the influence of the vector 
interaction and the neutral CSC phase , which favors the multiple critical-point structures or crossover for chiral 
restoration at low temperature. 

Of course, the produced u(d) quark vacuum constituent masses with different model parameters of the two-plus-one- 
flavor NJL model may range from 300-400 MeV, which are all phenomenologically acceptable just as the two-flavor 
case. Then, one can expect that all the critical-point structures found in Sec. |TT] should also appear in the two-plus- 
onc-flavor case, even considering the axial anomaly interaction term. 

IV. CONCLUSIONS AND OUTLOOK 

In this paper, we have explored the effect of the repulsive vector-vector interaction combined with electric-charge 
neutrality in [3 equilibrium on the chiral phase and CSC phase transitions within both two-flavor and two-plus-one- 
flavor NJL models. 

For the two-flavor case with the presence of the neutral CSC phase, we demonstrated that, with the help of the 
repulsive vector interaction, there always exists a parameter window in the NJL models which favors the appearance 
of a multiple chiral critical-point structure for a wide range of the vacuum quark mass, i.e., from 300 MeV to 400 
MeV. Besides the two- and three-critical-points structures found in 0, [2l[ and [2l[, respectively, we have shown 
for the first time that a four-critical-points structure of the QCD phase diagram can be realized ; such a multiple 
critical-point structure is caused by the joint effect of positive [i e and Gy ■ Because the dynamical strange quark mass 
is still relatively large near the boundary of the chiral transition, the multiple critical-point structures present in the 
two-flavor case also appear in the two-plus-one-flavor case. For the intermediate diquark coupling case, the number 
of critical points changes asl^2^4-^2^0 with an increasing vector coupling in the two-plus-one-flavor NJL 
model. In general, one can expect that different model parameters may possibly give other order of the number of 
the critical points as the vector coupling is increased. 

Although our analysis is based on a low-energy effective model which inherently has, more or less, a parameter de- 
pendence, we have seen that the physical mechanism to realize the multiple critical-point structure is solely dependent 
on the basic ingredients of the effective quark dynamics and thermodynamics. Therefore, we believe that the results 
obtained in the present work should be taken seriously and examined in other effective models of QCD, or hopefully 
lattice QCD simulations. Our result also has a meaningful implication for the study of phase transitions in condensed 
matter physics. That means some external constraints enforced on the system can lead to the formation or expansion 
of the coexisting phase, and the competition between two order parameters can give rise to multiple critical points. 
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Last but not least, we emphasize that we have shown for the first time that the repulsive vector interaction which 
should generically exist between the quarks does suppress the chromomagnetic instability related to the asymmetric 
homogeneous 2CSC phase. With increasing vector interaction, the unstable region associated with chromomagnetic 
instability shrinks towards lower temperatures and higher chemical potentials. This means that the vector interaction 
can partially or even totally resolve the chromomagnetic instability problem. 

Note that to cure the chromomagnetic instability, inhomogeneous asymmetric color superconductivity phases such 
as the LOFF phase and the gluonic phase were proposed in the literature. For the inhomogeneous phase, beside 
the condensate (ipjo'ip} , there is no reason to rule out the appearance of another new condensate , (ip"fip), when 
considering the vector interaction. The effect of both the timelike vector condensate and the spacelike condensate on 
the asymmetric inhomogeneous CSC phase will be reported in our future work[48|. 
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